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Abstract 

We consider non-perturbative QCD effects on the energy spectrum of either one of the 
photons in B — > X s ^^. These are due to the subprocesses in which a charm quark loop 
interacts with a self-consistently produced background static QCD field. The magnitude 
is estimated to be a few percents in B — * X s ^, but can be quite substantial in B s — > 77. 
An extension of the Euler-Heisenberg Lagrangian is given. 
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I. Introduction 

The flavor changing neutral current inclusive transition B — > X s ^ has been a subject of 
intense interest during the last few years. The basic theoretical framework is the Standard 
electroweak Model at scales ~ m w or m t . QCD short distance corrections [1] are then 
incorporated in it via renormalization group technique to yield an effective weak interaction 
Hamiltonian valid at scales ~ mj relevant for B decay processes. The special role that 
B — > Xs'j enjoys is related to its being rather clean and as such is much more model 
independent compared to exclusive B-decays. Its rate has recently been measured [2] with 
increased accuracy and is in remarkable agreement with theoretical estimates. This implies 
that revelation of New Physics will have to wait and that further confrontations have to be 
mounted. In light of these developments, it is natural to consider other inclusive channels 
which as a whole will separate out contributions from various operators in the effective 
Hamiltonian. The branching ratios may be somewhat lower in other novel processes, but 
they will still be amenable to experiments in the new facilities either under construction 
or being proposed for the not too distant future. One such process is the inclusive process 
B — > X s 77, which is expected to be ~ 10~ 2 smaller in its branching ratio relative to 
B — > X s -f. Just like B — > X s 7, however, it is relatively clean after some proper precaution 
to take out effects due to strong resonances such as r\ c at its peak to the two photon 
spectrum. It will then provide further opportunities for testing the whole technology of 
weak decays, or better yet in pointing towards some clues of New Physics. 

There have been some theoretical activities for the process B X s ^, which corre- 
sponds at the quark level to the transition b — > S77. Calculations were first done on the 
basis of pure electroweak theory [3-6] and subsequently improved to include the leading 
order renormalization group improved QCD effects [7-9]. These investigations (like most 
of the investigations for B — > X s ^) are mostly based on the free quark decay of b, the jus- 
tification of which is from the heavy quark effective theory (HQET) [10]. According to the 
argument given, corrections to the free quark results in inclusive processes are suppressed 
by powers of ( ) 2 . Recently, however, Voloshin [11] has shown that corrections which 
scale like ( ^ CP ) 2 should also exist. This last number is ~ 1 and thus has the potential of 
damaging all free quark estimates of inclusive B decays. For the leading process B — > X s ^, 
however, it has been shown that the Voloshin type of corrections have small coefficients 
multiplied to ( A ^ CP ) 2 and the overall effects are about 3% of the main term [12-15]. In 
fact, a more positive view towards this kind of corrections is to interprete them as system- 
atic accounts of the long distance non-resonant contributions of cc intermediate states. We 
shall subscribe to this constructive point of view and the present article is an investigation 
of the related effects in the parallel process B — > X577. We shall show that while the 
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corresponding corrections apparently have terms which scale like ( Q ^?' ) 2 ? where k is 
a typical photon momentum, in addition to corrections of the ( A ^ C£) ) 2 type, the overall 
effects are also only a few percents, just as in B — > X s ^. 

The plan of this article is as follows: in the next section, we shall briefly summarize 
some perturbative results of b — > S77 as a way to introduce our notation. In section III, 
we shall discuss the relevant diagrams which give rise to -K and -K corrections to this 
process. They come from a charm quark loop, from which an almost static gluon is emitted 
in addition to the two photons. Explicit formulaes will be given and matrix elements by 
HQET will be used to estimate their contributions to the decay amplitude. Some numerical 
work will be presented in the last section, followed by concluding remarks. 

II. b — > S77 in Standard Model with leading QCD corrections 

Radiative b — > s processes is best described in the framework of the following effective 
Hamiltonian 



where V^-'s are the CKM matix elements, Gf is the Fermi constant, Ci are the QCD 
improved Wilson coefficients, and the CVs are local operators: 
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Fig.l: One particle irreducible diagram (a) and one particle reducible diagrams (b) 
and (c) for b(p) — > s(p') n f(ki) r Y(k2)- Diagrams with (ki, fx) <-> (kz, v) should be added for 
(b) and (c). 



Upon using this effective Hamiltonian, the basic diagrams for the process b — > S77 are 
shown in fig.(l). Pure electroweak theory can be reclaimed by retaining only the C2O2 
term in fig. (la) and C7O7 term in figs.(lb-lc), with the given by Inami and Lim 
[16]. As expected, figs.(lb-lc) exhibit infra-red divergences, which are cancelled by the 
virtual radiative corrections to the process b — + S7 as far as any observables are concerned. 
Alternatively, a practical measurable rate for the b — > S77 can be formed by making a 
cut on the lower end of the energy of either of the two photons. This rate in the pure 
electroweak theory is known to be dominated by the one particle reducible diagrams of 
figs.(lb-lc); inclusion of QCD effects enhances their dominance, because C7 becomes even 
bigger. Therefore, in what follows we will assume that the principal diagrams for b — > S77 
are figs.(lb-lc), with only C7O7 insertions. We will compare the 'non-perturbative' gluon 
corrections to be discussed in the next section relative to this contribution. 

The b — > S77 amplitude generated by figs.(lb-lc) strictly speaking cannot be described 
by an effective local Hamiltonian. However, for a fixed value of the invariants p • ki, p-k^. 
and fci • &2, we can formally write an effective Hamiltonian of the form bT ^ u sA^ A u , where 

carries 7-matrices and momenta. After summing over the polarizations of the s-quark 
and the photons, the transition rate can be represented as < B\T{H\. Jl)H e f f(l))\B >, 
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where the following transition operator is defined 

e 2 G F X t Q d C 7 



< B\T(Hl ff (l)H eff (l)))\B >= ( -j§^Y < B\b(W}r^^(W 7 )^b\B >, 



(3) 



where X t = V t bV£ s and 



(W 7 )^ = - [- 2p.k 2 ^ ll ^ mhR + msL ^ mb + h)lv 



+ 



1 



-^ u {m s - pj + #i)#i7 M (m b .R + m s L) 



(4) 



2p' ■ k 2 
+ (ki, [i <-> k 2 , ^)]- 

A saliant feature of the spectrum generated by eqs.(3-4) is that as a function of the invariant 
mass of the two photons, it peaks at small values. This point will become useful in 
simplifying the one-gluon corrections, which we now immediately turn to. 

III. — , and —r corrections to b — > S77 

Because the coefficient C 2 is much bigger than 04,6, we shall discard effects due to 
O4 6 in our discussion. As for the operator 2 , we rewrite it by a Fierz transformation as 

1 



where 



Q 2 = -0i + 20i, 



0i = -c-i^—Lc-s^—Lb. 



(5) 
(6) 



Fig.la 





Fig. 2: (a) A charm quark loop which generates an effective bs^g vertex; (b) and (c) 
represent one particle reducible diagramss, which involve the bs^g vertex above. Diagrams 
with (ki, ffj <-> (k 2 , v) should be added. 



The operator 0\ generates a b-s transition with a c quark loop. Hooking onto it a 
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photon and a gluon, fig. (2a) leads to an effective bs^g Voloshin local effective vertex 



2 - Vol = 7^Jl P Y LbG ^ d " pPa ' (7) 

in the soft gluon limit, where F pa = ^e poiXK F\ K , e 0123 = i) i s the dual tensor to the 
electromagnetic field, and G a a ^ is the gluon field tensor. In this approximation which was 
used successfully and justified in ref. [11-15], the gluon inside the B-meson has been treated 
as a static field. The effective bs'jg- vertex in eq.(7) generates IPR b s'yyg amplitudes 
as shown in figs.(2b-2c). They cannot strictly be written as local effective interactions. 
However, as in the case of figs.(lb-lc) discussed before, for a given set of values for the 
momentum variables, the matrix element of figs.(2b-2c) (plus those with fci, \i <-> k<2.-,v) 
can be formally looked upon as arising out of a Hamiltonian 

^H-^K-^o,, (8) 



with (G af3 = g s G a a ^) 



R = sT^G^bA^AM, (9) 



and 



+ (fci, /U <-> k 2 , v). 

(10) 
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Fig. 3: A one particle irreducible diagram due to a charm quark loop, which generates 
an effective b — > 577(7 vertex. Five other diagrams due to permutations of^'s and g should 
be added. 



Finally, fig. (3) represents an irreducible diagram for b — > 577(7 transition generated 
by the charm loop, again via 0\ of eq.(6), which scales like ^4. Unlike the vertex in 
fig. (2a), only the vector part of 0\ now gives a non-vanishing contribution to it. Fig. (3) 
can be written as an effective Hamiltonian, which is basically the photon-photon scattering 
amplitude. The general expression for arbitrary fci, k<i and k$ is quite complicated [17], 
but, in the same spirit taken earlier, the gluon inside the B-meson can be treated as a 
background static field. We first form the gluon field tensor and then consider the limit 
with its momentum k^ — > 0. In this provision, we obtain 



x ( G s a « F^H^-f - JL) + (64^ - 
+ (1280^ + 672^ + 1)/^] 

40 TYl <z ~ 

+ G a K F K ^H^ u F v5 [— - 16-^Iqo 

+ (-512=£ - 224=| - |)/H 

qo 1 rn 2 A 

+ Gi"F SK F--H^ [( F + ^) + (-16-^ + T )/ 1 

+ (-384^ - 208^ - f )//„] 

11 1 m 2 

+ G a K Hfi K F^ u F^ u + 7^2^) - 

+ (-112^ - - »)/!„] 

1Q9 T7? 2 Ifi 

- (a A 9"-F"")(%F^)^Gf [— + (-96-^ + ¥ )I 1 m 
+ (-2304^- 1152^ 



where 

1 dx\dx 2 dxs8(l — x\ — x 2 — xs) 1, 9 . J km 2 + t + \/t. 

— = — In { — ) 

m 2 c + Xl x 2 t 2t ^Aml + t-y/i 

«S = 'Vy*. h„„ = 8„h; - 8„H°. (12) 

have been defined for t = 2fci • > 0. The prescription w? c — ie should be taken to 
analytically continue to the physical region t < 0. As we mentioned in the last section, the 
region of interest for the spectrum b — > S77 turns out to be when one of the momenta of 
the two photons becomes soft. Then the expression in eq.(ll) is further simplified into a 
local one 




+ Gj^F^ff^) + Gfi^F^F^)] , 

(13) 

which is in momentum space 

H eff (3) = {-i^^){-^^)KOm)llG 5K bA^k 1 )A^k 2 ), (14) 

where 

(0/fl)£ =^ x ^^((fcx + fc)V - (fei + h) V) + 14(^7^ - - ^2) 

+ 7^"(A;*^ - fcf - fc 2 ) • 7 + 3^((A; 1 + fc 2 ) V - (^ + fc 2 ) V)** • k 2 
- 7(g» 5 g™ - g^g^h ■ k 2 (h - k 2 ) ■ 7 
+ 3(0""^ - g^kl)k\k 2 • 7 + 3(^ 2 5 - g^k^h • 7 
+ S(g^k1 - g^kO^ki ■ k 2 + S(g vS k2 - g UK k d 2 )rh ■ k 2 
+ 7(g» s k% - g^k 6 2 )k\{k x + k 2 ) • 7 + 7(g v 'kZ - g^k^h + k 2 ) • 7 
+ U(g^k s 2 - g^k^h ■ k 2 + U(g v «k{ - g ud k^h ■ k 2 ]L. 

(15) 

The quarks inside the B-meson are in constant interaction with each other via soft 
gluon exchanges. The B — > X s 77 transition can thus be regarded as the quark transition 
b — > S77 in the presence of a background gluon field. This is the approach of Voloshin, 
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who further regarded this background gluon field as static. To first order in the strong 
QCD coupling, the effective Hamiltonians H e ff(2) of eq.(8) amd H e ff(3) of eq.(14), in 
which a zero momentum gluon is emitted from a c-quark loop, would then implement 
this dynamics. The corresponding amplitudes so generated would add coherently with the 
main amplitude of eq.(3). These correction terms are small (to be justified a posteriori), 
so that the addition to the principal transition operator would be the interference terms 
of amplitudes of fig. (2-3) with those of figs.(lb-c). 

The interference terms between H e ff(l) due to C 7 7 and H e ff{2) are 



^ RlTYw-t f,\TT / 9 ui Rn / e 2 Q d G F X t C 7 4GpAt e 2 Q d C 2 . 
< B\T(H eff (l)H eff (2))\B >=( 



(16) 



and the complex conjugate. Corresponding to the interference of H e ff(l) and H e ff(3), we 
have 



< B\T(H eff (l)H eff (3))\B >=( -^_)(^^-)(_^_) 

< B\Kwi)^^^(0 IR )ZG^b\B >, 

and its complex conjugate. Eqs.(3, 16 and 17) give the total transition amplitude 

< T >= < B\T(Hl ff (l)H eff (l))\B >+(< B\T(Hl ff (l)H eff (2))\B > +c.c.) 
+ ( < B\T(Hl ff (l)H eff (3))\B > +c.c.) 

We rely on HQFT to evaluate the matrix elements [18] 



(17) 



(18) 



< B(v)\bTb\B(v) >= lTr[l^r], (19) 

and 

< B(v)\bTG afi b\B(v) >= yTr[^ri^^], (20) 

where v is related to the momentum of b by p = m^v, V is any Dirac structure, and A2 is 

related to the B* — B mass splitting with a numerical value A2 = 0.12(Gev) 2 . We have 
normalized < B(v)\B(v) >= 1. The rate for b(p) — > s{p')^{ki)^{k2) is 

rp 1 X 4 ( / . , x d 3 p' d 3 h d 3 k 2 , . 

dT = 2W 5 {p ~ v ~ kl " k2) < T> • (21) 



IV. Numerical Results and Discussion 
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We would like to estimate the relative contributions of the three terms in eq.(18). 
As discussed earlier, we shall interprete corrections from the second and the third term 
as an indication of the long distance non-perturbative effects away from the peaks of the 
appropriate cc resonances. In this regard, the case for r\ c was addressed by the authors in 
ref.(8). They concluded that its effects on the spectrum of the two photon invariant mass 
are very localised. Our study here is on a different aspect of the same kind of issues. 

We find it convenient to present our numerical work in a quantity introduced by these 
authors. This is the spectrum of the photon with the lower energy and that of the photon 
with the higher energy, defined by 

^ - S ^k 9{±kMlk ' <22) 

where the integration domain, as in ref.(8), is restricted by the requirements that the 
energy of each photon be larger than E™ in = 100 Mev, and that the angles between any 
two outgoing particles be bigger than 20 degrees. These make it experimentally practical 
to distinguish from b — > sj. We fix \i = nib = 4.8Gev, m c = 1.5Gev and m s = 450Mei> , 
which lead to C 2 = 1.09 and C 7 = -0.31. Figs. (4-5) exhibit 50 x and 100 x the two 
interference terms separately vs. the main term of the spectrum from Eq.(18). Note that 
the photon energy k± is measured in unit of nib- Also, all the curves have been normalized 
to of the main term at k± = 0.086. The results show that despite the -^or^ scaling 

ClK\ " Tib TTb ^ 

the static gluon corrections from the charm loop is at most a few percents and hence not 
experimentally observable at this time. This is similar to the results for the one photon 
process, where however the corrections scale only as -K. 



Figure 4 Figure 5 
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Fig. 4-' The spectra of higher energy photon (HI) and lower energy photon 
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(LO) due to the first main term of Eq.(18) and 50x of the corresponding quantities due to 
the second term, denoted as HICOR and LOCOR, respectively, are shown. k\ is measured 
in unit of mi,. All curves have been devided by the value of LO at k\ = 0.086. 

Fig. 5: The spectra of higher energy photon (HI) and lower energy photon 
(LO) due to the first main term of Eq.(18) and lOOx of the corresponding quantities due 
to the third term, denoted as HICOR and LOCOR, respectively, are shown. k\ is measured 
in unit of nib. All curves have been devided by the value of LO at k\ = 0.086. 

The fact that a correction term which apparently goes like (^l 2 -) 2 , where k's are 
the typical photon momenta, is small is obviously because the two photon invariant mass 
stays small in the regions with appreciable rate. This situation would change dramatically 
for the exclusive process B s — > 77 where the invariant mass — 2k\ • /C2 = w 2 . One can 
expect contributions from fig. (3) to be of order (^-) 4 . On the other hand, B s — > 77 is an 
exclusive channel, and therefore the same technique applied here and in refs. [11-15] cannot 
be directly relied upon. We are presently engaged in evaluating these effects. 

Before closing, we would like to point out that eq.(13) is in agreement with the Euler- 
Heisenberg Lagrangian [19] and that eq.(ll) is an exact extension of the photon-phonton 
scattering amplitude to a situation when one of the photons (G^ in this case) is static 
while two on-shell photons can have any four momenta. The last photon is off-shell to give 
non-trivial kinematics. 
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